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ABSTRACT 

We have carried out a systematic search for close supermassive black hole binaries (with sub-parsec 
separations) among z < 0.7 quasars observed spectroscopically in the Sloan Digital Sky Survey. Such 
binaries are predicted by models of supermassive black hole and host galaxy co-evolution, therefore 
their census and population properties constitute an important test of these models. Our working 
hypothesis is that one of the two black holes accretes at a much higher rate than the other and 
carries with it the only broad-emission line region of the system, making the system analogous to 
a single-lined spectroscopic binary star. Accordingly, we used an automatic technique based on 
spectroscopic principal component analysis to search for broad H/3 emission lines that are displaced 
from the rest-frame of the quasar by \Av\ > 1,000 km s^ 1 (corresponding to binary periods and 
separations of P ~ few x 100 yr and a ~ few x 0.1 pc, respectively, for masses ~ 10 s M Q ). This 
method can also yield candidates for rapidly recoiling black holes since their spectroscopic signature 
is similar. Our search yielded 88 candidates, several of which were previously identified and discussed 
in the literature. The widths of the broad H/3 lines are typical among quasars but the shifts are 
extreme. We found a correlation between the peak offset and skewness of the broad H/3 profiles 
(there is an extended wing on the opposite side of the profile from the shifted peak), which suggests 
that the profiles we have selected share a common physical explanation. The general properties of 
the narrow emission lines are typical of quasars. We carried out followup spectroscopic observations 
of 68 objects to search for changes in the peak velocities of the H/3 lines (the time interval in the 
observer's frame between the original and new observations is 1-10 yr and 5.7-10.0 yr in 2/3 of 
the cases). We measured statistically significant changes (at 99% confidence) in 14 objects, with 
resulting accelerations between —120 and +120 km s _1 yr -1 . The above results, taken at face value, 
are broadly consistent with predictions for the number of close supermassive binaries in the Sloan 
Digital Sky Survey quasar sample. However, such a comparison is complicated by several theoretical 
and observational uncertainties, such as the fact that the observable we employ to select objects 
depends on a combination of several degenerate intrinsic parameters of a binary. We emphasize that 
interpretation of the offset broad emission lines as signatures of supermassive binaries is subject to 
many significant caveats. Many more followup observations over a long temporal baseline are needed 
to characterize the variability pattern of the broad lines and test that this pattern is indeed consistent 
with orbital motion. The possibility that some of the objects in this sample are rapidly recoiling 
black holes remains open as the available data do not provide strong constraints for this scenario. 

All tables and two large figures in this preprint have been abridged. A version 
with all the figures and a complete list of the sample objects can be found at: 

http://www2.astro.psu.edu/users/mce/preprints/SBHB.pdf (5 MB) 

Subject headings: galaxies: active — (galaxies:) quasars: emission lines — (galaxies:) quasars: general 
— line: profiles 



1. INTRODUCTION 
1.1. Theoretical Background and Motivation 
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The realization that the masses of the nuclear super- 
massive black holes (BHs) hosted by massive galaxies 
are related to the mass of the stellar spheroid of the 
host (e.g., iFerrarese fc Merrittl 120001: iGebhardt et al.l 
I2000T) has led to the development of new numerical 
and semi-anal ytic galaxy ev o lution models. Such 
mode ls (e.g.. IVolonteri et all 120031: IDi Matteo et al.l 
I2005t IHopkins et al.l 120061) attribute the observed 
correlation to co-evolution of the black hole and the 
host via hierarchical galaxy mergers and accretion. 
Close supermassive black hole binaries (SBHBs) ap- 
pear to be an inevitable stage in the evolution of 
the post-merger remnant. SBHBs have also been 
invoked to explain a variety of observations. Exam- 
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pies include stellar dynamical models that seek to 
explain the formation of cores in elliptical galaxies 
following; a merger and mass deficit s therein (e.g., 
IMilosavljevic et all [200l iMerrittl [2006h . and the mor- 
phologies of some radio sour ces (such as precessing jets 
and X- shaped source s; see IRoos. Kaastra. fc Hummel 
1993t iRomero et all l2000t iMemtt fe Ekersl 120021: 
Gopal-Krishna. Biermann. fe Wiital 12003ft . SBHBs are 
also expected to be among the primary sources of grav- 
itational waves detectable by the Laser Inter ferometer 
Space Antenna (LISA; e.g.. iSesana e t al. 200j|). 

The evolution of SBHBs following the merger 
of their parent galaxies was fi rst described by 
iBeeelman. Blandford. fc Reesl (|1980ft . The early evolu- 
tion (sinking of the BHs towards the center of the rem- 
nant and decay of the orbit via dynamical friction) pro- 
ceeds relatively quickly, on a time scale comparable to 
the dynamical time of the parent galaxies (~ 10 8 years). 
Once the binary "hardens" (i.e., the orbital speed be- 
comes comparable to the stellar velocity dispersion) the 
evolution slows down as the orbital decay relies primar- 
ily on scattering of stars. If the number of stars whose 
orbits allow close encounters with the binary can be re- 
plenished efficiently, the process can shrink the orbit to 
a separation where the emission of gravitational waves 
becomes an effective orbital decay mechanism and the 
binary can coalesce within a Hubble time. However, it 
is still unknown whether this replenishment is possible 
via stellar dynamical processes (see, for example, dis- 
cussi o ns in IMilosavljevic fc Merrittl [ 2003: Bc rczik et al.l 
[200l ISesana. Haardt. fc Madaul l2007t) . Therefore, the 
possibility that the orbital evolution of the binary may 
stall at a separation < 1 pc remains open (this is known 
as "the last parsec problem"). Alternative mechanisms 
for the decay of the binary orbit have been explored in 
an effort to solve this problem, most notably drag by 
a large gaseous reservoir (e.g. a disk) or torques on a 
binary at short separations by a circumbinary disk(e.g. 
Escala et all200llpotti et al.l2007l[2009bllCuadra et al.l 
2009HLodato et al.ll2009f ). 

Despite their importance and the theoretical atten- 
tion they have received, robust, dynamical evidence for 
the existence of large numbers of SBHBs at separations 
< lpc (the slowest ph ase in their evolution according to 
IBeeelman et al.lll980ft remains elusive. Discovering such 
systems and obtaining a census of their population prop- 
erties would serve as a test of galaxy evolution models 
and would provide valuable constraints for stellar and gas 
dynamical models for the decay of the binary orbit. It 
would also validate explanations of observed phenomena 
that appeal to such SBHBs. 

The above evolutionary scenarios along with the new 
observational work on this subject (described below) 
have prompted the development of models for the ob- 
servational signature of such binaries, such as the be- 
havior of the emission lines from gas accreting onto the 
two BHs. The observational signature that is thought to 
provide the most direct, dynamical evidence for SBHBs 
is based on an analogy with spectroscopic binary stars. 
In this picture the emission lines from gas bound to one 
or both of the BHs (the broad-line region or BLR) are 
Doppler-shif ted according to th e line- of-sight velocity of 
the source. IBogdanovic et all (|2008f ) have studied the 
interaction of the less massive (hereafter secondary) BH 



with the accretion disk around the more massive (pri- 
mary) one in a close, eccentric binary and have high- 
lighted the po tentially complex variability of t he em ission 
line profiles. lHavasak i, Mineshige. fc Sudul ([20071 ) and 
ICuadra et al.l (|2009D developed smoothed-particle hydro- 
dynamic (SPH) models for the gas flow from a circumbi- 
nary disk onto disks around the individual BHs, based 
on which they make predictions about the relative ac- 
cretion rates. In the context of such a picture of a bi - 
nary within a circumbinary disk, IMontuori et al.l ([20100 
have used photoionization models to compute the rel- 
ative strengths of the broad emission lines from these 
systems in a n effort to de v ise spe ctroscopic diagnostics 
for SBHBs. iShen fc Loebl (|2010aD present a toy model 
of the BLRs in SBHB in which they employ test parti- 
cles to represent clouds of gas in gravitationally bound 
orbits in the system and discuss its consequences for ob- 
servational searches for SBHBs. Although this model 
is unrealistic^ and the predicted line profiles cannot be 
used as a basis for quantitative tests, the authors make 
the very important point that, at short separations the 
BLR would envelope the entire binary. This point was 
al so made by M. Penston in 198 8 , as n oted in footnote 3 
of IChen. Halpern. fc FilipTjenkol (Il989h . 

A related development in numerical relativity is the 
calculation of the recoil speed (or "kick") of a BH re- 
sulting from the inspiral and merger of two BHs. A rc- 
vi ew of the recent li teratur e on this subject can be found 
in ICentrella et ahl (2010a. b). The most extreme recoil 
speed expected from this effect is ~ 4000 km s _1 for a 
special orientation of the spins of the two BHs relative to 
their orbital plane. However, under most other configu- 
rations, th e recoil speed is conside ra bly lower than thi s 
value (e.g- lCampanelli et al.ll2007allB:lBake"r et al.ll2008ft . 
As noted by IBogdanovic. Reynolds, fc Miller! (|2007|) , in 
a realistic astrophysical case the spins of the two BHs 
may be aligned with the orbital axis as a result of the 
previous evolution of the binary in a gaseous environ- 
ment. In such a case, the recoil speeds would be rather 
low, namely of o rder a few hundred km s . Indeed, 
iDotti et al.l (|2010f ) predict, based on detailed and exten- 
sive simulations, that the recoil speeds are likely to be 
distributed in the range 10-100 km s _1 . 

As in the SBHB case, finding recoiling BHs and study- 
ing their population properties would be quite valuable. 
It would provide a test of calculations for the recoil speed 
but also much-needed constraints on models of the as- 
trophysical implications for such kicks. A particularly 
interesting consequence of kicks is that they can pro- 
pel the BH out of the host galaxy if they are strong 
enough (i.e., > 1,000 km s _1 ), which could influence 
the c osmological evolution of BHs (e.g., iVolonteri et al.l 
2_010i). Otherwise, kicks can disturb the morphology 
of the host galaxy or off set an active nucleus from 
the center of its ho st (e.g.. IGualandris fc Merrittl [20081 : 
IKomossa fc Merrittl I2008t IBlecha et al.ll2007D . A poten- 
tial observational signature of recoiling BHs is a Dopplcr 

5 In addition to the shortcomings mentioned by the authors 
themselves, we also note that the basic premise of a BLR con- 
sisting of discrete clouds i s not supported by observations (see, for 
example, [Laor ct al. 2 0061 , and references therein). Even if a spher- 
ical system with crossing orbits like the one assumed could be set 
up, it would be destroyed by co llisions on a very short time scale 
(see IMathews fc Capriottill 19851) . 
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shift of the broad emission lines origi nating in a BLR 
that is bound to the recoiling BH (see lLoebl l2007D ■ An- 
other possibility is that the recoiling BH may accrete 
gas that it encounters along its trajectory, which could 
lead to a pro l onged AGN phase lasting up to 0.1-1 Gyr 
iBlecha et al.l (|2007f ). For large kicks, the spectroscopic 
signature of such a system is similar to the signature 
one might expect from some sub-parsec SBHBs, imply- 
ing that a systematic search for SBHBs based on this 
idea may also yield useful constraints on the frequency 
of large kicks. 

1.2. Observational Searches for Close Binary and 
Recoiling Holes in Historical Perspective 

Until recently, direct imaging had resolved very few 
(presumably unbound ) wide BH pairs, such as NGC 6240 



( Komossa ct a 



Man ess et al.l 



ll20"03h. and a close pair, CSO 040 2+379 



)03), 

200HRodriguez et al J 120061 120091 ). Re- 



cent spectroscopic surveys, however, have uncovered a 
substantial population of moderate-redshift AGNs with 
offset or double-peaked, narrow [O ill] emission lines. Al- 
though some of the offset or double-peaked [O in] lines 
coul d arise in biconical out flows from single active nuclei 
(see iCrenshaw et al.l I2010T ). some of them turn out to 
trace dual active nuclei in merging galaxies, as show n 
by followup st u dies (e.g.. i Comerford et al.l l2009alfbl ; 
Liu et al II 201 Oallbl : IShe^T et al.boiObHSmith et al.ll2010bl: 
Fu et, all f20Trf T These dual active nuclei have sepa- 
rations of order a few kiloparsec (a representative ex- 
ample is COSMOS J10 0043. 15+020637.2, reported by 
IComerford etld1l2009bn . 

In contrast, close, bound binaries, at sub-parsec sepa- 
rations, have been quite elusive as they cannot be re- 
solved spatially. Nevertheless, such close binaries are 
of particular interest because they represent the slow- 
est phase in the evolut ion of such a system, according to 
iBegelman et al.l l)1980D . In a well-known case, the ~ 12- 
year periodic modulation in the long-term light curve of 
OJ 287 has been expla ined in terms of a close SBHB (e.g., 
iValtonen et~aT1 120081 and refer ences therein). Other 
cases rely on a different type of indirect evidence, the 
detection of broad emission-lines with displaced peaks in 
quasar spectra. The underlying hypothesis here is that 
one or both of the two BHs in a SBHB has a gaseous 
BLR associated with it that emits the broad, permit- 
ted lines that are characteristic of quasar spectra. Thus, 
the orbital motion of the binary causes the lines to shift 
periodically making the spectrum analogous to that of 
a single- or double-lined spectroscopi c binary star (see 
iKomberd fl968l ; iBegj man et aLl I1980D . In the context 
of this hypothesis. iGaskelll (|1983l 1 19841) identified two 
quasars, 3C 227 and Mrk 668, whose broad H/3 emis- 
sion lines were shifted from their nominal positions by 
2,000-3,000 km s" 1 and proposed that these are exam- 
ples of close SBHBs in which only one of the two BLRs 
is visible. He also suggested that quasars with double- 
peaked emission lines, such as 3C 390.3, could repre- 
sent SBHB in which both BLRs are visible. Indeed, 
an observed systematic drift of the blue-shifted, broad 
H/3 peak in the spectrum of 3C 390.3 bet ween 1968 an d 
1988 was consistent with this hypothesis (Gaskcll 1996). 
However, this behavior is a rather uncommon feature 
of double-peaked emission lines (e.g.. iGezari et al.ll2007l : 
ILewis etall 12010ft and, in fact, the trend observed in 



3C 390.3 did not continue p ast 1990 (|Eracleous et al.l 
Il997t iShapovalova et al.ll2001[) . Moreover, in a detailed 
study of the long-term var iability of double-peak ed emis- 
sion lines in three quasars. lEracleous et al.l (|1997ft did not 
detect the behavior expected from a binary BLR associ- 
ated with an SBHB and also noted a number of con- 
tradictions between this interpretation and the observed 
properties of these objects. Finally, long-term variabil- 
ity studies of about two dozen doub le-peaked emitters 
(jGezari et al.ll2007t ILewis et all 12010ft do not show the 
signature of orbital motion and a number of arguments 
suggest that broad double-peaked lines are much more 
li kely to originate in an acc retion disk (see discussion 
in lEracleous fc Harpernll2003ft . Nevertheless, the binary 
black hole hypothesis remains a reasonable explanation 
for broad, single-peaked Bahncr lines that arc shifted 
from their nominal wavelengths. It is interesting that 
a number of such cases have been known for quite some 
time among nearby S eyfert galaxie s. An inspection of the 
spectroscopic atlas of lStirpel (|1990ft reveals several exam- 
ples, of which Mrk 876 and NGC 55480 arc particularly 
striking. 

The recent theoretical work on the evolution of black 
holes in merging galaxies has been accompanied by new 
observational work directed at identifying close (sub- 
parsec) SBHBs. A handful of quasars from the Sloan Dig- 
ital Sky Survey (SDSS), whose broad H/3 lines have peaks 
that are displaced from their nominal wavelengths by a 
few thousand km s , have been noted and discussed. 
The quasar SDSS J092712.65+294344.0 was interpreted 



as a SBH B bv iBog danovic. Eracleous. fc Sigurdssonl 



(|2009ft an d iDotti et al.l (l2009aft. al thoug h it was pro 



posed by iKomossa. Zhou, fc Liul (|2008f ) as a candi- 
date recoiling BH (see also our discussion below). 
Th e quasar SDSS J 1 53636 22+044127.0 was identified 
by iBoroson fc Lauerl ([2009) who attributed its strong, 
broad, blucshiftcd H/3 peak to a BLR around one of the 
me mbers of a SBH B. However, iChornock et al.l (|2010l ) 
and IGaskelll (|2010l ) called this interpretation into ques- 
tion based on the presence of a prominent red shoulder in 
the Ha profile, which is also discernible in the H/3 profile, 
and suggested that the Balmcr lines originate in a per- 
turbed disk around a single black hole. Yet another inter- 
pretat ion for this object was offered bv lTang fc Grindlavl 
(|2009ft . who proposed that this quasar does harbor a 
SBHB and the profiles of the Balmer lines comprise con- 
tributions from an accretion disk around the primary 
BH, a BLR associated with the secondary BH, and a 
more extended gaseous region surrou nding the binary. 
An in terpretation similar to that of ITang fc Grindlavl 
(|2009T ) was also discussed by IBarrows et alJ (|2011l ) for 
the H/3 profile of SDSS J093201.60+031858.7. A fourth 
case, 4C+22.25 (a.k.a. SDSS J100021. 80+223318.7) 
was r ecently reported and discussed by iDecarli et al.l 
(|2010ft . The H/3 profile of this object is reminiscent of 
SDSS J092712.65+294344.0 (a broad peak blueshifted by 
8,700 km s" 1 ), but the large blucshift of the peak makes 
a rapidly-recoiling BH interpretation untenable. In a 

6 The variations of the H/3 profile of NGC 5548 were discussed 
by IPeterson. Korista. fc Cotal I|1987I 1 in the context of a SBHB 
model. However, the wealth of data available today indicate that 
the variability pr operties of this object do not support such an in- 
terpretation (e.g. , Wanders & Peterson 1996; Peterso n et al.ir2002l ; 
ISergeev et~aT]|2007t l . 
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very recent study, iTsalamantza et aLI (|2011f) have used 
an automatic method to search systematically through 
the SDSS DR7 spectroscopic database (including 54,586 
objects classified as quasars and 3,929 objects classified 
as galaxies at 0.1 < z < 1.5) for objects with displaced 
emission line peaks. They have uncovered five new SBHB 
candidates, in addition to the ones noted above. All of 
the above objects should be regarded as preliminary can- 
didates for SBHBs in the absence of additional evidence. 
This is because our lack of understanding of the structure 
and properties of the BLR does not preclude that such 
unusual line profiles can originate from the BLR around 
a single BH. The case would be strengthened if the dis- 
placed peaks are observed to drift in time, as one would 
expect for the lines of a spectroscopic binary star. 

A systematic search for recoiling BHs was undertaken 
only recently, mot ivated by the n e w dev elopments in nu- 
merical relativity. iBonning et al.l (|2007| ) studied the dis- 
tribution of shifts of the broad H/3 lines in the SDSS 
DR3 quasar sample and found that shifts greater than 
800 km s _1 occur in only 0.2% of cases, with progres- 
sively more stringent limits for larger shift s. The lat- 
est st udies for the distribution of kicks (e.g.. iDotti et al.l 
120101 ) suggest that the vast majority of recoil speeds 
should be <100kms _1 , which means that recoiling 
BHs would not be e asily identifiable in the sample of 
IBonning et all (|2007f ). 

Three candidate rapidly recoiling black holes have 
been noted and discussed in the recent literature, but 
their cases are ambiguous or controversial. As we men- 
tioned above, SDSS J 92712 65+294344.0 was identi- 
fied by iKomossa et al.1 (j2008f ) via its shifted H/3 line, 
who suggested that this is a rapidly-recoiling black hole. 
However, a number of different interpretations were of- 
fered by followup papers , including a binary black hole 
(|Bogdanovic et al.l 120091: IDotti et al.l l2009af) . a chance 
supe rposition of objects at different redshifts on the 
sky ([Shields et al.1 l2009bl ) , and a super position of ob- 
jects in the same, massive galaxy cluster (jHcckma n et al.l 
l2009h . Similarly SDSS J105041.35+345631.3 was discov- 
ered through a visual inspect ion of a large numb er of 
SPSS quasars and d iscussed bv lShields et al.l (12009a! ) and 
iSmith et al.l (|2010a[ ). These authors noted that the rel- 
atively symmetric displaced peaks of the Balmer lines 
suggest a rapidly recoiling BH but differences between 
the profiles of broad lines from different ions (H I and 
Mg n) suggests that the lines are actually double-peaked 
but extremely asymmetric, which would p oint to an ori- 
gin in a perturbed accretion disk. Finally, iCivano et al.l 
(|2010( ) suggested that one of the two nuclei CID-42 (a.k.a. 
COSMOS J100043. 15+020637.2; see the first paragraph 
of this section) could also be a rapidly recoiling BH. 

The ambiguity surrounding the objects mentioned in 
the previous paragraph stems from the fact that the 
available data do not allow us to definitively distinguish 
between SBHBs and rapidly recoiling BHs. The detec- 
tion of a change in the velocity offset of the broad lines 
would favor the SBHB over the recoiling BH interpreta- 
tion, although it would not provide iron-clad evidence, as 
we discuss in more detail in later sections of this paper. 
One can argue against a recoiling BH interpretation in 
some particular cases on theoretical grounds, since recoil 
speeds far in excess of 100 km s _1 are extremely unlikely 



(e.g., IDotti et al.ll2010D and speeds above 4,000 km s" 1 
are unattainable (e.g., iCampanelli et al.l l2007allb |) . An 
indirect test has b een described by Bonning et al.l (|2007T ) 
and re-iterated bv lShields et al.l ()2009af ) in which the dis- 
placed broad lines have fairly symmetric profiles and the 
prominent narrow lines have similar widths, regardless of 
ionization state. The underlying hypothesis here is that 
the broad lines are emitted from an otherwi se ordinary 
BLR that is attached to the recoiling BH (see lLoebll2007j) 
and that the narrow lines are emitted from the galaxy 
where the BH originated, which is now illuminated from 
the outside. 

2. SCOPE AND GOALS OF OUR PROGRAM AND 
CONSIDERATIONS UNDERLYING THE SELECTION OF 
CANDIDATES AND SEARCH STRATEGY 

Motivated by the above, we have embarked on a sys- 
tematic search for quasars whose broad (single-peaked) 
H/3 lines are offset from their nominal wavelengths by 
\Av\ > 1,000 km s _1 . We regard such objects as candi- 
dates for close (sub-parsec) SBHBs or rapidly recoiling 
BHs and have begun followup observations to search for 
changes in the shifts of their broad lines. Even though 
we adopt this as our working hypothesis, we emphasize 
that such line profiles need not be signposts of SBHBs. 
One can envision scenarios in which such line profiles can 
originate from the BLR associated with a single BH that 
does not have a companion. Thus, the mere detection 
of an offset in the broad Balmer lines yields only weak 
candidates for SBHBs or rapidly recoiling BHs. The de- 
tection of gradual drifts in the broad Balmer lines via 
followup observations strengthens the case for a SBHB, 
but a convincing demonstration that a specific object 
harbors a SBHB requires us to observe (ideally) a few 
orbital cycles. We return to these issues in SJ7J where we 
offer a more extensive discus sion. 

Based on the models of iVolonteri et al.l (|2009D we 
would optimistically expect up to ~ 160 sub-parsec SB- 
HBs in the sample that we have searched. However, 
there are significant theoretical uncertainties and models 
by these authors based on different assumptions yield a 
number that is lower by an order of magnitude. Since the 
SBHB phase is thought to be considerably longer than 
the accretion phase of a recoiling BH and since the ex- 
pected distribution of recoil speeds is heavily weighted 
towards low speeds (\v\ <C 1, 000 km s _1 ) we expect that 
only a very small fraction of the t he quasars we select will 
be candidate reco i ling B Hs (see iBogdanovic et al.l 120071 
IDotti et al.lf2009al f2010h . 

If the binaries evolve within a large-scale gaseous 
disk, as suggested by recent theoretical studies (e.g., 
IDotti et al.ll2007l I2009bf ). we would expect that their or- 
bits will have circularized by the time they reach sub- 
parsec orbital separations. At this evolutionary stage, 
binaries open a hole at the center of the circumbinary 
disk and the secondary BH orbits closer to the gas reser- 
voir and has easier access to it. Therefore, the accretion 
rate onto the secondary is considerably higher than that 
onto the primary (see the recent work by Hayasaki et 
al. 2007 and Cuadra et al. 2009, building on previous 
calculations by Artymowicz & Lubow 1996 and Gould & 
Rix 2000), which leads us to attribute any shifted emis- 
sion lines we detect to a BLR associated with the sec- 
ondary. The resonant interaction of the binary with the 
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circumbinary disk could lead to an increase in the ec- 
centricity. This is a somewhat uncertain process, which 
depends on the thermodynamic properties of the cir- 
cumbina ry disk and the amount of gas accreted by the 
BHs ( seelArtvmowicz et aHll99lHArmitae;e fc Natarajanl 
120051 iCuadra et all 120091 iRodig et all 120111 ). We con- 
sider the accreting secondary scenario to be more likely 
on theoretical grounds, but our technique does not pre- 
clude us from finding systems i n which the BLR is as- 
sociated with the primary (e.g. iChang et alJ l2010h . or 
systems with a BLR that envelops the entire binary. 

Since the main observable we employ is the velocity 
offset of the broad H/3 line, we consider here how the 
properties of detectable binaries might depend on this 
quantity. For the sake of s implicity we assume that the 
orbit is circular. Following Bogdano vic et all ()2009f ). we 
write the magnitude of the observed velocity shift (i.e., 
projected along the line of sight) of the broad lines from 
the BLR of the secondary, u 2 , in terms of the true orbital 
speed, V 2 , as u 2 = V 2 sini | sin0|, where i is the inclina- 
tion of the orbit (i = 0° is face on), (j> = 2ir(t — to)/P 
is the orbital phase, and P is the period. We express 
the period and orbital separation of the binary in terms 
of the total mass, M 8 = (M x + M 2 )/10 8 M Q , mass ra- 
tio, q = M2/M1 < 1, and the projected velocity of the 
secondary, u 2> 3 = u 2 /10 3 km s _1 as: 



and 



P = 2652 M 8 



0.432 M 8 



sin 1 sin ( 



(1 + q) m 2 ,3 



sin i sin i 



(1 + q) u 2 ,3 



yr 



pc. 



(1) 



(2) 



If the emission lines trace the motion of the the primary 
instead, the above equations can be adapted by replac- 
ing U2,3 by U\fi/q, where 1*1,3 = iti/10 3 km s _1 {u\ is 
the projected speed of the primary). For the sake of es- 
timating typical binary periods and separations, we set 
i = (j> = 45°, which yields 



P 



332 M a 



(1 + q) 3 \ sin 45° sin45 c 



and 



0.11 M« 



(1 + 9) 2 



x 2,3 



sini I sin0 
sin 45° sin 45° 



yr 



pc. 



(3) 



(4) 



In other words, if we measure the projected speed of the 
secondary BH to be of order 10 3 km s _1 in a binary with 
5 « 1 and Mg ~ 1, then it is reasonable to expect the 
binary to have a period of order a few centuries and a 
sub-parsec separation. We note here that the mass of 
the secondary cannot be arbitrarily small since it is con- 
strained by the requirement that the accretion luminosity 
be high enough for this object to appear as a quasar. If 
the emission lines trace the motion of the primary in- 
stead, then the period and separation would be smaller 
by factors of q 3 and q 2 , respectively. These long orbital 
periods highlight a significant difficulty in securely iden- 
tifying SBHBs: orbital periods are very long compared to 
a human lifetime, hence the observation of even a single 
orbital cycle requires considerable patience and persis- 
tence. In an optimistic scenario, where the binary has a 



mass and a projected velocity offset such that Mg w 0.1 
and U2,3 ~ 2, the orbital period could be as short as a 
decade, however. We will return to this issue in $7] where 
we will consider the likelihood of observing a system with 
a specific set of properties. 

As we noted above, the case for a displaced broad line 
signaling the presence of a binary becomes stronger if 
the shift of the line is observed to change with time. In 
the context of the binary picture outlined in the previous 
paragraph, the projected acceleration is given by 



dt 



= 2.4 



4 3 (i + g) 3 

M 8 sin 3 i 



km s 1 yr 1 



(5) 
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"2,3(1 + l) 3 f sin 45 



M 8 

cos 01 / sin 45 



cos 45° V sin 



sin i 

4 

km s _1 yr^ 1 



The time intervals (in the observer's frame) between the 
original SDSS observations and the new observations car- 
ried out in our program are in the range of 5-10 years. 
Thus, for sub-parsec binaries with orbital periods of or- 
der centuries, it is reasonable to expect the line shifts 
to have changed by 100-200 km s _1 . Our survey was, 
therefore, designed to detect changes of this order. If, 
however, the broad line shifts trace the motion of the 
primary, then, for a given observed velocity the accel- 
erations are larger by a factor of q~ 3 compared to the 
expression in equation ([3]) and are easier to detect. 

It is noteworthy that \diij/dt\ cx Vj|cos</>|, while 
Uj oc Vj I sin 0| (where j = 1,2 labels the BH), i.e., the 
speed and acceleration are j cycle out of phase during 
their oscillation. This means that binaries near conjunc- 
tion should display small velocity offsets but large accel- 
erations, while binaries near quadrature should display 
lar ge velocity offsets and sm all accelerations. As noted 
by iBogdanovic et all (|2009f) . this represents an inherent 
difficulty in making a case for a SBHB because large ac- 
celerations should be accompanied by small line shifts 
and vice versa. 

With these considerations in mind, in fj3]we describe in 
detail how we selected our sample out of the SDSS spec- 
troscopic database and we present the basic properties 
of the resulting objects, while in i|4]we present a detailed 
study of the spectroscopic properties of these objects. In 
we present our followup ( "2nd-epoch" ) observations 
and compare the followup spectra to the SDSS spectra. 
The techniques we used to search for velocity changes in 
the displaced broad peaks are described in S}|5] and in Sj7] 
we summarize and discuss our results. 

3. TARGET SELECTION 

The basis for our identification of candidates for the 
target list is a principal co mponent analysis (PCA) tech- 
nique described in detail by |Boroioa^Tauei] ([20l3). We 
applied this technique to the SDSS spectra of 15,900 
QSOs having r e dshift z < 0.7, from the catalog of 
ISchneider et all (|2010D . We note that this catalog in- 
cludes only objects more luminous than Mj = —22. The 
input variables were the spectra themselves, shifted to 
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rest and trimmed to the wavelength range 4290-5400 A, 
centered roughly on the H/3 region. Eigenspectra were 
generated from a subset of these (1,000 objects) with the 
highest signal-to-noisc. Although the spectrum of almost 
every object could be modeled to within its noise by using 
a large number of eigenspectra, our goal was to identify 
the objects with the worst fits as candidates for our sam- 
ple. Thus, we reconstructed each of the 15,900 spectra 
using only the five most significant eigenspectra. This 
effectively limits the reconstruction such that it can only 
reproduce the most common features and characteristics 
of the line profiles. 

We adopted a threshold of x 2 /^ > 3 (where v is 
the number of degrees of freedom, i.e., the number of 
spectral pixels) as a starting point for identifying spec- 
tra that were fitted poorly and were therefore candidate 
offset-peak objects. This gave us 910 objects. Note that 
this threshold value is arbitrary and is dependent on the 
signal-to-noise ratio (S/N) of the spectrum. A very high- 
S/N spectrum with a slightly unusual line profile will 
result in a poor fit. However, we thought that this ap- 
proach would tend to err on the side of including many 
objects that we would later reject, and so it is inherently 
conservative. 

We visually inspected the 910 spectra identified by the 
automatic selection algorithm and applied the following 
additional criteria. We excluded most objects in which 
the peak of the broad H/3 line was shifted by less than 
1,000 km s _1 from the peak of the narrow H/3 line 
(interpreted as the systemic redshift) since any offset 
broad-line peak is blended with the narrow H/3 line 
and difficult to discern. Identifying the peak of the 
broad line was subjective in some cases, particularly in 
noisy spectra. We excluded H/3 profiles with shelves 
and inflections, if they had a dominant broad peak 
at v w (examples of quasars with such profiles 
that were excluded are SDSS J121716.08+080942.0, 
SDSS J143452.45+483942.8, SDSS 

J154348.62+401 324.9, and SPSS J 154929.43+023701.1; 
see Figure 5 of iZamfir et al.l l2010T l . We also excluded 
clear double-peaked emitters, i.e., objects with two 
clearly separated broad peaks on either side of the 
narrow H/3 line. This is because variability studies have 
tested and rejected t he hypothesis that such profiles are 
indicative of SBHBs (lEracleous et ai]|1997t iGezari et al.l 
[20071 iLewFeTaTl I2010D ." However, we retained objects 
whose broad H/3 profiles have a strong displaced peak 
plus an extended wing or shoulder on the opposite 
side of the narrow H/3 line (a prime example of such 
an object is SPSS J15 3636.22+044127.0, presented by 
iBoroson fc Laueil I2009L while additional examples from 
our collection are identified in ^4.3|F 1. We also retained 
flat-topped profiles, if these were not centered at v = 
since they could be a manifes tation of a single bro ad-line 
region bound to a SBHB (see IShen fc Loebll2010afl . This 

7 Objects with such profiles have been termed "extreme double- 
peaked emitters" and have been discussed as altern ative interpre- 
tation s for broad displaced peaks to SBHBs (e. g. [Shields et al. 
I2009bl ; iLaiier fc Borosonll2009l ; IDecarli et aLll2010r) . These two ex- 
planations are not mutually exclusive, of course. The unusual 
profiles of extreme double-peaked emitters could be the result 
of perturbation of the line-emitting dis k from a companion BH 
in a SBHB (e.g lEracleous et all jl995l ; IBogdanovic et al.l [20081 ; 
I'l'ang &; Grin dlav 200E])! However, similar perturbations can be 
caused by the self gravity of the disk. 



subjective, visual classification process resulted in a 
sample of 88 objects that met all our cri teria. To these 
we ad ded two more objects listed by iStrateva et al.l 
(|2003|) . which had not been selected by our x 2 / v > 3 
condition. Again, we stress that this is not a complete 
sample in any sense, though we believe that it includes 
almost all of th e objects that fi t our c riteria. 

According to iVolonteri et al.l (|2009f ) there may be up 
to 150 detectable SBHBs with q > 0.01 in the sample 
of 15,900 quasars we have searched. However, models 
based on different assumptions predict up to an order 
of magnitude lower numbers, and the predictions are, in 
fact, consistent with zero when theoretical uncertainties 
are considered. 

The 88 objects making up the resulting sample, along 
with their basic properties are listed in Table [1] This 
sample is considerably larger than the few objects with 
case studies published so far. In columns 1-5 of this ta- 
ble we give the redshift, apparent V magnitude, Galactic 
V-band extinction in the direction of the source, and the 
absolute V magnitude for each quasar. These quanti- 
ties were obtained as described in i j4.ll Hereafter, we 
will adopt abbreviated designations for these objects us- 
ing the first six digits of the right ascension (i.e., down 
to integer seconds, e.g., Jhhmmss). With this conven- 
tion we can refer to the objects in Table [T] uniquely. In 
columns 6-10 of Table [TJ we list the particulars of the 
spectroscopic observations of each target, starting with 
the original SPSS observations (additional rows refer to 
followup observations described in f|5]). 

Some of the objects in Table [T] were identified previ- 
ously and discussed by other authors, as follows. 

• The quasar J140700 is Mrk 668, one of the firs t 
two candidate SBHBs proposed by iGaskelll (|1983ft . 
It is noteworth y that the other candidate SBHB of 
IGaskelll (p83h . 3C 227, was also observed during 
the SPSS (designated SPSS J094745. 14+072520.5) 
but was not selected by our algorithm. This is 
because the profile of its broad H/3 line at the 
time of the SPSS spectroscopic observation (UT 
2003/03/26) was relatively symmetric and centered 
at the nominal wavelength of the line. We will re- 
turn to these two objects in our discussion of inter- 
esting cases in later sections of the paper. 

• We have recovered three objects that have re- 
ceived attention in the recent literature: J092712 
and J105041 (see JL2]), as well J153636 (origi- 
nally discovered by this m ethod). We have als o 
recovered an object noted bv lBonning eEaD([2Q03), 
J0918333. In gj2] we compare in more detail the 
properties of the broad H/3 line s from our sample 
to the findings of iBonning et al.l ()2007T ) and to sta- 
tistical studies by other authors. 

• Our sample ha s 13 objects in common w ith the list 
of 32 objects of lTsalamantza et al.l (|2011[ ). who also 
used an automatic selection algorithm to search 
through the SPSS PR7 spectroscopic database. Of 
these 13 objects, four are the objects listed in the 
previous paragraph, two are newly selected SBHB 

8 IBon ning ct al. (2003) drew attention to this object because it 
had the largest velocity offset in their sample 
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TABLE 1 

Log of Spectroscopic Observations 

















Exposure 




Rest-Frame 


Object Name 




my 
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My 


Observation 


Instr. 


Time 




Wavelength 


SDSS J 




(mag) 


(mag) 


(mag) 


Date (UT) 


Code 6 


(s) 


S/N c 


Range (A) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 
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3756 


29 


3093-7505 
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002444.10 + 003221.3 


0.4024 
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-24.86 


2000/12/22 


S 
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2709-6594 
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015530.01 - 085704.0 


0.1648 


16.84 


0.080 


-22.65 


2001/09/16 


S 


2702 


32 


3263-7917 












2009/12/17 


M 


3600 


55 


3375-6524 


020011.53 - 093126.2 


0.3602 


17.96 


0.081 


-23.46 


2001/08/28 


S 


8104 


34 


2795-6780 












2009/12/17 


M 


3600 


28 


2890-5586 


021259.59 - 003029.5 


0.3945 


17.87 


0.117 


-23.83 


2000/09/29 


S 


2700 


29 


2714-6613 












2009/12/16 


M 


3600 


13 


2803-5433 


022014.57 - 072859.1 


0.2137 


18.65 


0.078 


-21.47 


2001/09/10 


S 


4501 


11 


3133-7598 












2009/12/16 


M 


3600 


16 


3220-6241 



The remaining objects have been omitted 
The full table can be found in the complete preprint 



a The rcdshift of the object as measured in this paper from the peak wavelenth of [O in] A5007 line. 

k The apparent V magnitude, determined from the SDSS PSF magnitudes, as described in fj^of the text. 

c The Galactic visual extinction, taken from Schlcgcl, Finkbcincr, & Davis (1998) 

The absolute V magnitude, computed as described in JHof the text. 
c The telescope and instrument configuration, as described in Table [2] 

^ The signal-to-noise ratio (S/N) in the continuum near the line of interest. When the spectrum includes the H/3 line we 
give the S/N in the continuum near this line, at 4600 A. If the spectrum includes only the Mg II A2S00 line, we give the 
S/N in the continuum near this line, at 2900 A. 



TABLE 2 

List of Telescopes and Instruments 



Instrument 






Rcsolution a 


Code 


Observatory, Telescope, and Spectrograph 


Spectral Elements 


(A) 


S 


SDSS 2.5m, SDSS spectrograph 


640, 440 mm -1 grisms, d = 3" fibers 


2.7 


M 


MDM, Hiltner 2.4m, Boiler k, Chivens CCD spectrograph 


G350L grating (150 mm^ 1 ), l'.'O slit 


7.6 


K 


KPNO, Mayall 4m, Ritchie-Cretien spectrograph 


KPC-007 grating (600 mm" 1 ), l'.'B slit 


3.1 


Pb 


Palomar, Hale 5m, Double Spectrograph 


blue arm: 600 mm -1 grating, 1"5 slit 


1.0 


Pr 




red arm: 600 mm -1 grating, 1'.'5 slit 


3.1 


HI 


Hobby-Ebcrly Telescope, Low-Resolution Spectrograph 


Gl grism (300 mm" 1 ), l'/5 slit 


14.5 


112 




G2 grism (600 mm- 1 ), 1'.'5 slit 


5.6 



candidates by both us and them (J115449 and 
J171448), while the remaining seven are objects 
that we regard as SBHB candidates while they 
have classified as having asymmetric line profiles 
or "ot her." There are fou r SBH B candidates in the 
list of lTsalamantza et al.l (j 2 1 If ) that are not in our 
sample: J0932+0318, J1000+2233, J1012+2613, 
and J1539+3333 (using their naming convention). 
The second and third of these four they also classi- 
fied as double-peaked emitters (we agree with this 
classification) which is why they are not included 
in our sample. 

• Five of the objects in our sample, J001224, 



J015530, J094603, Jl 1 4 755, and J120924 were 
noted by iZamfir et al.l (|2010l ) as unusual cases 
that illustrate the wide diversity of Balmer line 
profiles in quasars. Several other objects in the 
IZamfir et al.1 (j2010[ ) study are worthy of mention, 
because they are not included in our sample. Two 
objects, J101912 and J 16 2345, are below t he lu - 
minosity threshold of the Schneid er et al.l (|2010f ) 
QSO list. Three objects, J074948, J084203, and 
J143511, have broad H/3 lines with peaks that are 
less than 1,000 km s -1 from the narrow H/3 line. 
One of these appeared in our candidate list with 
a x 2 l v value of 3.92 and was rejected upon visual 
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inspection. The other two had x 2 /^ values just 
below our cutoff. 

• Nine of the objects in our sample are also in the 
sample o f z < 0.33 quas a rs wit h double-peaked Ha 
lines by iStrateva et al.l (|2003l ). selected from the 
SDSS DR3. These are J001224, J022014, J074157, 
J081329, J093844, J110742, J115644, J143455, and 
J172711. Of these, J081329 and J172711 have 
Balmcr-linc profiles whose centroid is somewhat 
blueshifted and with a strong blue peak and a 
weak red shoulder. As such their classification is 
somewhat ambiguous. The remaining seven ob- 
jects have very skewed Balmer-line profiles with 
one clearly displaced peak, but a second peak is 
not clearly present, making their classification as 
double-peaked emitters questionable. Thus, we 
have retained all of these objects in our sample. We 
note that J001224 was use d as an illustra t ive exa m- 
plc of SBHB candidate bv lShen fc Loebl (|2010al) . 

In the next section we present the basic properties of 
the sample objects and we study their spectra, paying 
particular attention to the profiles of their H/3 lines. 

4. PROPERTIES OF THE SAMPLE 

4.1. SDSS Observations and Basic Properties of Sample 

Objects 

The apparent V magnitudes of the sample quasars 
were determined from the SDSS PSF magnitudes in the 
q and r bands, us ing the transformation equations of 
IJester et al.l (|2005i ). We obtained the redshifts of the 
sample quasars by measuring the peak wavelength of the 
[O ill] A5007 line and by adopting a rest vacuum wave- 
length of 5008.239 A for this line (since the SDSS spectra 
are on a vacuum wavelength scale). From this redshift, 
we dete rmined the d istance modulus using the prescrip- 
tions of Hogg ( 1999 adopting the following cosmolog- 
ical parameters: Hq = 73 km s _1 Mpc -1 , Qm = 0.27, 
and = 0.73. Thus, we obtained the absolute V mag- 
nitudes, taking into account Galac tic extinction (from 
iSchlegel. Finkbeiner. fc Davisl 119981) . but not K correc- 
tions. The redshift, apparent V magnitude, and absolute 
V magnitude distributions are shown in the histograms 
of Figure [T] We note that the median redshift of our 
targets is 0.32, while the full range of redshifts extends 
from 0.077 to 0.713. 

In columns 6-10 of TablcQ]we gi ye a log of the spectro- 
scopic observations of the objects in our sample, includ- 
ing the UT date of the observation, the telescope and 
instrument configuration, the exposure time, the signal- 
to-noise ratio in the spectrum (S/N) , and the rest- frame 
wavelength range covered by the spectrum. The table 
includes a separate block for each object, with the first 
row of each block giving the particulars of the SDSS 
observations. The instrument configurations are sum- 
marized in Table [2] It is noteworthy that the spectral 

9 We used a program written by Ben- 
jamin Weiner, which is publically available at 

http://mingus.as.arizona.edu/~bjw/software/ We com- 

pared the luminosity distances trom this program to 
those of Ned Wright's cosmology calculator, available at 
http : //www. astro .ucla. edu/7,7Ewright/CosmoCalc .html, and 
found them to agree within 0.1%. 



resolution of the SDSS spectra varies with wavelength 
such that the velocity resolution improves linearly with 
the log of the wavelength (see the discu ssion and illus- 
tratio n in Appendix B and Figure 16 of iBernardi et al~l 
|2003| ). Thus, in Table [2] we list the spectral resolution 
at 6400 A, which corresponds to the location of the H/3 
line for the median redshift of our sample. The S/N was 
determined from the fluctuations about the mean contin- 
uum level in a 50 A wide window near the emission line 
of interest. In particular, when the spectrum includes 
the H/3 line we give the S/N in the continuum near this 
line, at 4600 A (Av « -16,600 km s" 1 relative to H/3). 
If the spectrum includes only the Mg n A2800 line, we 
give the S/N in the continuum near this line, at 2900 A 
(Av « +10, 500 km s^ 1 relative to Mg n). 

The spectra of the Mg II, H/3, and Ha lines of the 
objects in our sample are shown in Figure |5|on a com- 
mon velocity scale. The broad-line spectra shown in this 
figure illustrate that it is easier to find broad displaced 
peaks in the H/3 profiles, especially when the displace- 
ment is small, because the peak of the broad Ha profile 
is "contaminated" by the narrow Ha+[N n] complex. For 
example, in J130534 and J140007, one can easily discern 
an offset broad peak in the H/3 but not in the Ha pro- 
files. Vice versa, the Ha profiles show more clearly the 
shape of the red side of the profile where the H/3 line may 
suffer from severe contamination by the [O in] doublet. 
For example, in J075403, J091928, and J121113, what 
appears to be an extended red wing in the H/3 profile ac- 
tually has the form of a shoulder or weak red peak in the 
Ha profile The profiles of the Mg II lines arc very similar 
to the profiles of the Balmer lines. Specifically, there is 
no discernible misalignment of the peaks of the different 
lines with only one exception: J105203. One must keep 
in mind, however, that the contamination of the peak of 
the broad line by the narrow doublet is more severe than 
in the case of H/3 and comparable in severity to the case 
of Ha. Moreover, since the Mg n lines are resonance 
lines, the profiles of the broad emission lines are often 
contaminated by associated absorption lines, which can 
be fairly strong. 

4.2. Processing of SDSS Spectra and Emission- Line 
Measurements 

To quantify the spectroscopic properties of our targets, 
we measured the properties of some of the strong narrow 
emission lines as well as the properties of the profiles of 
the broad H/3 lines from the SDSS spectra. 

Specifically, we measured the integrated fluxes and full 
widths at half maximum (FWHM) of the following nar- 
row lines: [Ne v] A3426, [O n] A3726,3729, [Ne m] A3869, 
H/3 and [O III] A5007. The widths and relative strengths 
of the narrow lines, especially the high-ionization ones, 
can serve as tests of the rapidly recoiling BH h y pothe - 
sis, as discussed, for exam ple, by iBonning et al.l (|2007lK 
iKomossa et~aT1 (|2008f ) and lShields et al.l (|2009air To iso- 
late the lines, we fitted a low order polynomial to the 
local effective continuum. In the case of the H/3 and 
[O ill] A5007, we used a continuum model that also in- 
cludes the contribution of the Fe II lines, as we describe 
later in this section. The line fluxes were then deter- 
mined by integrating the observed line profile, while the 
widths were determined by fitting a Gaussian to the line 
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Fig. 1. — Distribution of the basic properties of the objects in our sample. The redshifts were determined from the wavelengths of the 
[O in] A5007 line. The apparent V magnitudes correspond to the PSF magnitude obtai ned from SDSS images. The absolute V magnitudes 
take into account Galactic extinction towards the source. Further details are given in £|4,ll of the text. 



profile and then correcting the FWHM for the finite reso- 
lution of the spectrograph at the observe d wavelength of 
the li ne (see Appendix B and Figure 16 of lBernardi et al.l 
I2003D . We made no effort to de-blend the two lines in 
the [O n] AA3726,3729 doublet, therefore we do not re- 
port their widths. However, we did measure the width of 
the blended doublet, which allows us to discern in some 
cases whether the [O n] lines are narrower than other 
forbidden lines. The observed line fluxes were corrected 
for Galactic extinction using the values of A y listed 
in Ta ble [T] and assuming the extinction law of iSeatonl 
(|1979|) . The resulting measurements are reported in Ta- 
ble [3] where we give the line fluxes relative to that of the 
[O in] A5007 line before and after extinction corrections, 
the observed and corrected [O in] A5007 flux, and the 
corrected [O ill] A5007 luminosity, as well as the FWHM 
of the lines. 

In order to isolate the profiles of the broad H/3 lines we 
first subtracted a model of the underlying continuum and 
Fe II multiplcts, which was fitted to the spectral range 
4,000-5700 A. The underlying continuum consisted of 
a linear combination of a starlight template and a fea- 
tureless power law while the Fe II multiplets were mod- 
eled by suitably broaden ing a template, as described in 
iBoroson fe Greenl (|1992ft . As a last step we subtracted 
the narrow H/3 and [O in] lines by using the profile of the 
[O in] A5007 line as a template for the other two lines. As 
a cross-check, we used the profile of the [O in] A4959 as a 
template to subtract the [O III] A5007 line. The resulting 
broad H/3 profiles are shown on a common velocity scale 
in Figure [3] 

After isolating the broad H/3 line profiles we quantified 
their shapes by measuring some of their central moments 
and related quantities. The n th central moment of a line 
profile, [in, is defined in terms of the wavelengths and 
flux densities of the discrete pixels in the line profiles (A^ 
and fi respectively) as fi n = KJ2 i^i ~ W)" /«' wrierc 
K is a normalization constant defined by 1 / K = Y] fi 
and (A) = K ^2 A,; fi is the first moment or centroid of 
the line profile. Using the second and third central mo- 
ments we determined the standard deviation and skew- 
ness of the line profiles as a = fa' 2 and s = fa > 
respectively. We also determined the Pearson skewness 
coefficient, defined as p = ((A) — A Tn )/er, where A m is 
the median wavelength^- In Table [4] we report the cen- 
troid velocity shift relative to the nominal wavelength of 
H/3 using the relativistic Doppler formula, the velocity 

10 Because of the definitions we adopted for s and p, s oc —p. 



dispersion of the line, computed from the standard de- 
viation as ccr / Ao (where Ao is the nominal wavelength of 
H/3 and c is the speed of light), the skewness, and the 
Pearson skewness coefficient. 

In addition to the above quantities, we also measured 
the velocity offset of the peak of the broad line and its 
FWHM, which we also report in Table @] In order to 
make these measurements, we had to locate the peak of 
the broad line, which is somewhat subjective since the 
observed profiles have a finite S/N. Thus we attempted 
to determine the location of the peak of the broad line 
profile by fitting it with a Gaussian. This procedure leads 
to an uncertainty because the symmetric Gaussian does 
not always fit the line peak well, especially when the 
peak is asymmetric. To quantify this uncertainty, we 
made several measurements of the same line, varying the 
region around the peak used for the fit. In the end we 
adopted the maximum and minimum peak wavelengths 
from our trials to define the range of possible values and 
took their average as the best estimate of the peak wave- 
length. Thus, in Table [4] we report the velocity offset 
of the peak, computed from the relativistic Doppler for- 
mula, as well as its uncertainty. Finally, the FWHM was 
determined based on the height of the peak and expressed 
as a velocity width in the same manner. 

4.3. Demographics and Statistical Properties of the 
Broad H/3 Profiles and the Narrow Lines 

Figure @] illustrates the distribution of shifts of the 
peaks of the broad H/3 lines among objects in our sam- 
ple. The dearth of objects in the range ±500 km s _1 is 
a result of the fact that we purposely excluded most ob- 
jects with such small shifts. In the same figure we show 
the range of shifts (±3 standard deviations about a mean 
of 100 km s _1 ) measured by iBonning et all (|2007f ). who 
studied the H/3 profiles of 2598 quasars at 0.1 < z < 0.81 
drawn from a sample of approximately 13,000 objects 
from the SDSS DR5 and found 9 cases of large offsets, 
from « 1,000 km s _1 to a maximum of 2,667 km s . 
Our study appears to complement and supplement theirs 
since we recover a substantially larger number of highly 
offset peaks: 32 objects with \Av\ > 1500 km s _1 , a fac- 
tor of 8 larger than what we would predict just by scaling 
the sample sizes. Similar con clusions follow if we com- 
pare our results with those of IZamfir et a l. (2010), who 
studied the high-S/N spectra of 469 quasars with z < 0.7 
from the SDSS DR5. The vast majority of our objects be- 
long to their populati on B (FWHMha > 4 , 000 km s" 1 , 
weak Fe II lines; see iSulentic et al.l 120021) . whose dis- 
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Fig. 2. — ABRIDGED Mg II, H/3, and Ha spectra of our target objects. For each object we show a 3-pancl figure with the profiles of 
all three lines on a common velocity scale. The original SDSS spectra are plotted in black, while the post-SDSS spectra are color-coded as 
follows: KPNO 4m spectra are in red, Palomar 5m spectra are in blue, MDM 2.4m spectra are in green, and HET spectra are are in orange. 
For the purpose of this illustration, some spectra were binned to increase the S/N per bin. The vertical scale was adjusted specifically 
for the original SDSS spectra. Post-SDSS spectra have b een scaled and shifted so that the profiles of the broad lines match those of the 
SDSS spectra as closely as possible; see equation (|6j and £|5.2l of the text. As a result of this scaling, there is an apparent difference in the 
strengths of the narrow lines, if the flux of the broad line has varied between epochs. 
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TABLE 3 

Properties of Narrow Emission Lines 



Obs. Flux Relative to [O m] a Corr. Flux Relative to [O m] b FWHM (km s" 1 ^ 
Obs. Corr. 



Object 


[Ne v] 


[O II] 


[Nc m] 


H/3 


-F[oni] a 


[Ne v] 


[O ii] 


[Ne III] 


H/3 


T? b 
-HO III] 


-k[Oin] c 


[Ne v] 


[Ne III] 


H/3 


[O m] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


J001224 


<0.055 


0.462 


0.346 


0.128 


3.59 


<0.058 


0.483 


0.359 


0.129 


3.95 


5.73 




910 


510 


540 


J002444 


0.200 


0.130 


0.242 


0.193 


7.77 


0.207 


0.133 


0.247 


0.194 


8.19 


44.0 


890 


910 


420 


420 


J015530 


0.157 


0.166 


0.184 


0.237 


16.9 


0.162 


0.171 


0.189 


0.238 


18.0 


12.5 


670 


550 


470 


290 


J020011 


<0.105 


0.140 


0.220 


0.215 


3.58 


<0.109 


0.143 


0.225 


0.215 


3.77 


15.5 




1150 


720 


800 


J021259 


0.054 


0.132 


0.086 


0.075 


31.5 


0.056 


0.137 


0.089 


0.075 


34.0 


174. 


430 


490 


150 


440 


J022014 


<0.112 


0.564 


0.151 


0.219 


3.90 


<0.116 


0.579 


0.155 


0.220 


4.14 


5.14 




560 


400 


370 



The remaining objects have been omitted 



a Integrated line fluxes are expressed relative to the integrated flux of the [O III] A5007 line in units of 10 — 15 erg cm -2 s~ 1 , as measured from the 
observed spectrum, without any corrections. The uncertainty is 5% or less, unless otherwoise noted. 

k Integrated line fluxes are expressed relative to the integrated flux of the [O III] A5007 line in units of 10 _15 erg cm -2 s _1 , after correcting for 
Galactic extinction, as described in fl4.2| of the text. The uncertainty is 5% or less, unless otherwoise noted. 

c The luminosity of the [O III] A5007 line in units of 10 37 erg s — 1 , after correcting for Galactic extinction (see *|4.2| of the text). 

^ The full width at half maximum of the emission lines, corrected for the finite resolution of the spectrograph as described in EI4.2I of the text. 

The value is omitted when the line is not detected. 

c The uncertainty in these measurements is 10-15% 

* The spectrum of J092712 includes two sets of narrow lines therefore we report te properties of each set in separate rows. The first row gives the 
properties of the set that is blucshiftcd and appears to be at the same redshift as the broad lines. The second row gives the properties of the set 
that is rcshiftcd and has no broad lines associated with it. 



TABLE 4 
Properties of Broad H/3 Profiles 









Peak 


Centroid 


Velocity 




Pearson 




FWHM 


FWQM 


Velocity Shift 


Velocity Shift 


Dispersion 


Skewness 


Skewness 


Object 


(km s~ - 1 ) 


(km s -1 ) 


(km s" 1 ) 


(km s" 1 ) 


(km s —1 ) 


Coefficient 


Coefficient 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


J001224 


3590 


8630 


-1870 ± 80 


280 


3530 


0.253 


-0.197 


J002444 


9320 


13450 


-720 ± 150 


-80 


3720 


0.113 


-0.056 


J015530 


7320 


9180 


1430 ± 130 


610 


3100 


-0.039 


0.045 


J020011 


7660 


12490 


1650 ± 60 


460 


3570 


0.333 


0.025 


J021259 


7990 


11570 


-2300 ± 40 


-10 


3830 


0.323 


-0.106 


J022014 


8210 


10560 


-470 ± 90 


950 


3240 


0.495 


-0.102 



The remaining objects have been omitted 



a The redshift of the object as measured in this paper from the peak wavclenth of [O III] A5007 line. 

The apparent V magnitude, determined from the SDSS PSF magnitudes, as described in j4]of the text. 
c The Galactic visual extinction, taken from Schlcgcl, Finkbcincr. & Davis ( 1998) 

The absolute V magnitude, computed as described in j4]of the text. 
e The telescope and instrument configuration, as described in Table [2] 

^ The signal-to-noise ratio (S/N) in the continuum near the line of interest. When the spectrum includes the H/3 line 
we give the S/N in the continuum near this line, at 4600 A. If the spectrum includes only the Mg II A2800 line, we 
give the S/N in the continuum near this line, at 2900 A. 
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Fig. 3.— ABRIDGED Profiles of the broad H/9 lines from the 
SDSS spectra, shown on a common velocity scale. The narrow H/3 
and [O in] lines and underlying continuum were subtracted as de- 
scribed in 8^4.21 of the text. The resulting profiles were normalized 
arbitrarily and offset vertically from each other for clarity. The ver- 
tical dashed line shows the location of the narrow H/3 line, while the 
vertical dotted lines identify a window of ±5, 000 km s _1 from this 
line. The noise often seen at +8875 km s — 1 is a result of imperfect 
subtraction of the [O III] A5007 line. The He II A4686 line, which 
was not subtracted, is sometimes discernible at —10,800 km s — 1 . 



tribution of broad H/3 shifts shows some preference for 
rcdshifts and has extrema of ±3100 km s _1 . Thus, the 
objects we have identified here occupy the wings of the 
distribution of shifts of the peak of the broad H/3 line. 

Perusal of Figures [5] and [3] leads us to classify the 
Balmer line profiles into three broad families accord- 
ing to their shapes. The first family includes profiles 
whose widths are smaller than average and visual in- 
spection does not show them to be particularly asym- 
metric. Two of the objects discussed extensively in the 
recent literature, J092712 and J105041, arc prime ex- 
amples of objects in this family. Additional examples 
include J082930, J095036, J110050, J122811, 132704, 
and J162914. The second family includes skewed pro- 
files showing an extended wing or a shoulder on the side 
opposite from the direction of their displacement (e.g., 
profiles with a blue-shifted peak with an extended red 
wing or a weak red shoulder). Representative examples 
of this family with blue-shifted peaks (and red asymme- 
tries) are J001224, J091928, J115449, J123001, J125809, 
J133432, J152942, J153636, while examples with red- 
shifted peaks (and blue asymmetries) include J015530, 
J093653, J094603, J113651, J130534, J143455. The third 
family includes objects with fiat-topped profiles, such as 
J020011, J031715, J074157, J095539, J112751, J154340, 



Oh 

0) 
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-5000 5000 

Peak Velocity Shift (km s" 1 ) 

Fig. 4. — Distribution of the velo city shifts of the peaks of the 
broad H/3, measured as described in i|4.2l of the text. The dearth of 
objects in the rage ±500 km s _1 is a selection effect resulting from 
our inability to detect such small shifts with our method. The 
horizontal error b ar at the top of the frame indicates the range 
of shifts found by Bonning et al. (2007): ±3 standard deviations 
about a mean of 100 km s — 1 . 

J163020. 

In addition to the visual classification of line profiles 
we have also looked for trends in their quantitative prop- 
erties and we found a correlation between their shift and 
their skewness. We illustrate this correlation in Figurc[5j 
where we plot the Pearson skewness coefficient vs the 
peak shift. A similar correlation is recovered when we 
consider other pairs of quantities that track skewness and 
shift, such as the skewness coefficient and the centroid 
shift. As the shift of the lines increases, the profile be- 
comes more skewed. Moreover, the asymmetry manifests 
itself as an extended wing or weak shoulder on the side 
of the line profile that is opposite from the direction of 
the shift, with the result that the sign of the skewness 
tracks the sign of the shift. This correlation quantifies 
the trend we noted in the second family of line profiles 
in the previous paragraph. We have not found any cor- 
relation between the FWHM of the lines and the shift 
of their peaks nor between any other pairs of qu antities 
that track width and shift. iBonning et al.1 (|2007h do re- 
port that H/3 profiles with large shifts also have large 
FWHM but the vast majority of objects in their sam- 
ple show considerably smaller H/3 shifts than the objects 
in our sample and the shifts they regard as large are 
considerably smaller than the shifts we measure. There- 
fore, our results are not at odds with theirs. We did 
find a weak correlation between the width and the skew- 
n ess. These correlatio ns are consistent with the findings 
of lZamfir et al.l (|2010ft . shown in their Figures 6a and 6c. 

In Figures |5] and [7] we present the distribution of rel- 
ative intensities and FWHM of the the narrow lines us- 
ing the measurements presented in Table 03 The rela- 
tive intensity distributions resemble very closely those 
of Seyfert 1 an d Seyfe rt 1.5 galaxies in the compilation 
of iNagao et al.l (|2001[ ). which suggests that the ioniza- 
tion parameter and the shape of the ionizing contin- 
uum are similar to those of Seyfert galaxies. The dis- 
tribution of FWHM of the narrow line s is similar to 
what is observed in Seyfert g alaxies (see iWhittlel 119851 : 
iMoore. Cohen. &: Marcvlll996 ) and is also eviden t in th e 
composite quasar spectrum of lVanden Berk et al.l (|200lD . 
The average FWHM of the [O in] line is very similar to 
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Fig. 5. — Relation between the skewness of the broad H/3 profiles 
and the shift of their peaks. The Pearson skewness co effic ient and 
the shift of the peak were measured as described in 114.21 and are 
tabulated in Table |I] 



that found by iSalviander et "aD (|2007| ) in their study of 
the spectra of 1736 quasars from SDSS DR3. Both the 
[Ne v] and [Ne in] lines in our sample are broader than 
[O in], on average, and the [Ne v] lines are always as 
broad as or broader than the [O in] line within errors. 
This trend follows the correlatio n between FWHM and 
critical density note d for Syeferts (jFilippenko &: Halpernl 
119841: 1 Whittld fl985h . Thus the objects in our collection 
do not appear to stand out from typical active galactic 
nuclei as far as the properties of the narrow lines are 
concerned. 

5. FOLLOWUP OBSERVATIONS AND NEW SPECTRA 
5.1. Data Acquisition and Reduction 

We carried out followup observations of 68 objects from 
our sample using four different telescopes: the brightest 
objects were observed with the 2.4m Hiltncr telescope 
at MDM observatory, objects of intermediate brightness 
were observed with the 4m Mayall telescope at Kitt 
Peak National Observatory and the 5m Hale telescope 
at Palomar Observatory, while the faintest objects were 
observed with the 9.2m Hobby- Eberly Telescope (HET). 
With the followup observations we specifically targeted 
the H/3 lines of the target objects, although we also did 
attempt to obtain spectra of the Ha and Mg II lines 
whenever possible. The spectrographs and configura- 
tions we used for the observations are summarized in 
Table [21 along with the spectral resolution that each con- 
figuration can attain. The log of observations is included 
in Table [T] where, in addition to the observation date and 
exposure time, we also report the rest-frame wavelength 
range covered by each spectrum and the S/N achieved 
in the continuum near the line of interest (measured as 
described in fj3j, above). The weather conditions were 
varied during the followup observations and some spec- 
tra were taken through clouds. As a result, there are 
instances where we obtained a lower S/N than the origi- 
nal SDSS spectrum of the same objects in a comparable 
exposure time. 

The spectra were reduced in a standard manner using 



IRAI0 for the first steps of the reduction process (bias 
subtraction, flat field division, sky subtraction, and ex- 
traction of raw spectra) and our own programs for the 
subsequent calibration (wavelength calibration, correc- 
tion for continuous atmospheric extinction and discrete 
absorption bands, and flux calibration). Accurate wave- 
length calibration is important for our purposes since 
we aim to look for shifts in the H/3 lines between the 
SDSS and post-SDSS spectra; thus we give here the de- 
tails of the process. To obtain the transformation be- 
tween pixel number and observed (air) wavelength we 
fitted a polynomial to the table of arc line pixel loca- 
tion vs wavelength. Depending on the instrument, we 
used 20-65 arc lines and achieved a fit with root-mean- 
squarc (r.m.s.) residuals of better than 0.10 pixels with a 
polynomial of order 4 or lower. Thus, the relative wave- 
length scale is good to 10 km s -1 for the KPNO 4m and 
Palomar 5m spectra, 13 km s" 1 for the HET spectra 
and 18 km s _1 for the MDM 2.4m spectra. In compar- 
ison, the relative wavelength scale of the SDSS spectra 
is good to 0.07 pixels or 5 km s _1 (r.m.s.). After the 
spectra were fully calibrated, we applied heliocentric cor- 
rections (amounting to 0.15 A or less) and converted the 
wavelength scale from air wavelengths to vacuum wave- 
lengths, following the convention for SDSS spectra. For 
the last step we used t he relation betw een vacuum and air 
wavelengths given by iMortonl ()1991[ ) , which yields cor- 
rections between approximately 0.5 and 2.5 A over the 
wavelength range of our spectra. With these corrections 
the r.m.s. shift between the SDSS and post-SDSS spec- 
tra amounted to 0.67 A or 41 km s _1 , determined by 
cross-correlation of the spectra over the region around 
the [O in] AA4959,5007 lines. Thus, in a final step we re- 
fined the alignment of spectra by rectifying these residual 
shifts. The final alignment is limited by our ability to de- 
termine the shift of the [O in] lines by cross-correlation. 
In about 80% of the cases the r.m.s. alignment uncer- 
tainty is 7 km s _1 or better, while in the other 20% of 
cases it lies between 7 and 25 km s _1 . 

5.2. Qualitative Comparison of New and Old Spectra 

The distribution of time intervals (in the observer's 
frame) between the original SDSS observations and the 
followup observations is shown in Figure [5] The median 
of the time interval distribution is 7.0 years, with 70% of 
the objects having a time interval between observations 
of at least 5.7 years. 

In Figure [2] we overplot the post-SDSS spectra on the 
original SDSS spectra for comparison. Since our primary 
goal is to look for changes in the profiles of the broad 
emission lines, we have scaled and shifted the post-SDSS 
spectra so as to match these profiles as closely as possible. 
We applied a transformation to the flux density scale of 
the form 

fx = afx + b x , (6) 

where a and b\ arc constants, determined by minimizing 
the differences between the SDSS and post-SDSS spec- 
tra in selected regions of the broad line profiles and the 

11 IRAF, the Interactive Reduction and Analysis Facility, is dis- 
tributed by the National Optical Astronomy Observatory, which 
is operated by the Association of Universities for Research in As- 
tronomy (AURA) under coo pera tive agreement with the National 
Science Foundation (see xlOhttp:/ /iraf.noao.cduT). 
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Fig. 6. — Distribution of four different narrow-line ratios among t he q uasars in our sample, after correction for Galactic extinction. The 
measurements were made from the SDSS spectra, as described in £14.21 of the text, and the results are tabulated in Table \3\ The shaded 
portions of bins represent actual measurements while the hollow portions represent limits. 



adjacent continuum. This scaling convention highlights 
changes in the profiles of the two spectra but it hides 
changes in the flux of the broad line and the continuum. 
Changes in the flux of the broad line can still be dis- 
cerned by examining the intensities of the narrow lines 
(most notably [O in] ), which are not expected to change 
in the time interval between observations. In other words 
an apparent mismatch in the strength of the narrow lines 
between two spectra indicates that flux of the broad line 
has changed between observations. 

An inspection of Figure [2] reveals changes in the broad 
H/3 lines between the SDSS and post-SDSS observations 
in many cases. Changes in the profile shapes are the easi- 
est to discern. Examples of objects displaying this type of 
variability include J093653, J112751, J113706, J140700, 
and J143123. There are also several cases of changes in 
the integrated flux of the line, discernible through mis- 
matches in the strengths of the [O in] lines, including 
but not limited to J021259, J094603, J111329, J115449, 
J132704, and J143455. Yet another category of changes 
comprises cases where the broad line has shifted in wave- 
length/velocity by a small amount. The most obvious 
example of such a change is J095036, while J094603, 
J140251, and J180545 represent examples of more subtle 
shifts of this type. This type of variability if, of course, 
of particular interest here, therefore we study it further 
in 36l below. 



SEARCH FOR VELOCITY CHANGES IN THE BROAD H/3 
PROFILES 



6.1. 



A x 2 Cross-Correlation Method for Determining 
Velocity Changes 

To determine any changes in the velocities the offset 
peaks of the broad H/3 lines between the SDSS and post- 
SDSS spectra we developed an algorithm that is a variant 
of the standard cross-correlation method. In our method 
we shift one of the two spectra in small steps and at each 
step we compare it with the other spectrum via the x 2 
test. Based on several tests, we prefer this method over 
the standard cross-correlation (which employs the over- 
lap integral rather than the % 2 test at each step) for the 
following reasons: (a) it allows us to check whether the 
profile of the H/3 line has changed between epochs (this is 
important because small changes in the profile can mimic 
a velocity change) , (b) the x 2 minimum is much sharper 
and deeper than the overlap integral maximum, thus eas- 
ier to locate, especially if we are relying on the wings of 
the broad line for our comparison, (c) it allows us to de- 
termine uncertainties and limits on the shift in a more 
straightforward and less computationally intensive way, 
without the need for extensive simulations for each pair 
of spectra being compared. 

Our x 2 cross-correlation method is applied to a pair 
of spectra as follows: before any comparison is made, a 
linear transformation is applied to the flux scale of the 
second spectrum of the pair following equation ([5]) , in or- 
der to match the continuum and the broad H/3 line. This 
transformation creates a mismatch in the narrow lines, of 
course, but it allows us to carry out a direct comparison 
of the broad lines, without subtracting the continuum or 
the narrow lines, which would introduce additional sys- 
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Fig. 7. — Distribution of the FWHM of high-excitation narrow 
lines of quasars in our sample, after correction for the finite resolu- 
tion of the spectrograph. The mea surements were made from the 
SDSS spectra, as described in i|4,2l of the text, and the results are 
tabulated in Table \3\ The upper panel shows the distribution of 
widths of the [O III] A5007 lines, while the lower panel shows the 
distributions of [Ne III] A3869 (black, solid line with hollow bins) 
and [Ne v] A3426 (shaded bins with no outline). 
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Fig. 8. — Distribution of the time intervals (in the observer's 
frame) between the original SDSS spectra and the followup spectra 
obtained here. The median of this distribution is 7.0 years. For 
70% of the objects the time interval between observations is at 
least 5.7 years. 

tematic errors. The wavelength scale of the second spec- 
trum is then shifted in small steps (typically 0.1-0.2 A), 
rcbinned to the wavelength scale of the first spectrum, 
and finally the x 2 statistic is evaluated to compare the 
profiles of the broad H/3 lines in one or two windows that 
are not affected by narrow lines. These windows are cho- 
sen to encompass the offset broad peak and/or the steep 
wings or shoulders of the line so as to maximize the sen- 
sitivity to small shifts. The final output is a x 2 curve as a 
function of shift. The minimum of this curve, determined 




by fitting a parabola to the five lowest points, represents 
the optimal wavelength shift. The 99% (or "2.6 cr") con- 
fidence interval about the minimum (for one interesting 
parameter) is determined by finding the shifts that corre- 
spond to Xmin + 6.63 (see lLampton. Margon. fc Bowverl 
[19761) . In Figure 1 we show graphically the application of 
this technique in three example shift that is sta- 

tistically significant, a shift that is not statistically signif- 
icant, and an apparent shift that is caused by a variation 
in the broad H/3 profile between the two epochs. 

We verified through simulations that xLin + 6-63 does 
indeed reflect the 99% confidence interval. Each simula- 
tion was carried out by creating a copy of an observed 
high-S/N spectrum, adding synthetic noise, and then 
measuring the shift between the two spectra. We re- 
peated this test for several examples of observed spectra 
from our collection. In 1000 realizations in each case 
we found that about 1% of the optimal shifts departed 
from zero by more than the amount predicted by the 
Xmin + 6.63 condition, as expected. An additional check 
of the procedure was carried out by switching the order 
of the spectra supplied to the \ 2 cross-correlation algo- 
rithm and comparing the resulting shifts. Indeed the 
distribution of discrepancies between the two shifts, nor- 
malized by the 99% confidence error bar was found to 
be approximately Gaussian with a standard deviation of 
1/2.6, as expected. 

6.2. Application of the x 2 Cross- Correlation Method 
and Results 

We applied the above method to identify variable ve- 
locity candidates which we then scrutinized carefully to 
ensure that the shifts were genuine. We rejected cases 
where substantial variations of the broad line profiles 
mimic shifts (see example in Figure ■ We screened the 
initial candidates as follows: First we checked that the 
result was the same regardless of which of the two spectra 
was shifted and rebinned in the cross-correlation proce- 
dure (if more than one followup observation was carried 
out we also compared the results from the different obser- 
vations). Then we inspected the spectra to ensure that 
the shift could be identified visually. In cases where the 
shift was subtle, we experimented by varying the spec- 
tral windows over which the x 2 statistic was computed 
in order to verify that the result was robust. Finally, for 
the candidates that passed the above tests, we carried 
out simulations to verify the uncertainties: we made a 
copy of the spectrum with the higher-S/N added syn- 
thetic noise to match that of the \owei-S/N spectrum 
and then applied the cross-correlation algorithm. After 
1000 realizations we examined the distribution of shifts 
about zero to verify that the statistical uncertainty in 
the shift was not underestimated. After arriving at the 
final statistical uncertainty in the velocity change of each 
object, we added it in quadrature with the uncertainty in 
aligning the SDSS and post-SDSS spectra via the [O in] 
lines (see the last paragraph of §5.1[) to obtain the total 
uncertainty. 

In Table [5] we summarize the results of theta cross- 
correlation analysis. In 14 cases the broad H/3 peak ve- 
locities change significantly (at 99% confidence). In 38 
cases we obtain only limits to possible shifts (at 99% 
confidence). In 10 of these cases we encountered small 
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Fig. 9. — Examples of the application of our cross-correlation technique (see i|6.1l of the text) to three pairs of spectra. In the lower panel 
of each set we show the the two spectra appropriately scaled and superposed for comparison. The original SDSS spectrum is shown in black, 
the post-SDSS spectrum (which is the one that is scaled, shifted and rcbinncd by the cross-correlation algorithm) is shown in red, while 
the residuals after subtraction of the two is shown in blue (in this illustration the residuals represent the difference between the un-shifted 
spectra). Since the scaling of the two spectra is such that the continuum and broad lines match, there is typically a mismatch between 
the narrow lines, which indicates that the continuum and/or the broad lines have varied between the two epochs. The vertical dotted lines 
mark the nominal positions of the narrow H/3 and [O III] lines. The horizontal error bars show the sections of the broad line profiles used 
by the cross-correlation algorithm. The upper panel of each set shows the variation of x 2 with the shift applied to the post-SDSS (red) 
spectrum. The red arc at the trough of the \ 2 curve shows the global minimum and the (blue) horizontal line shows the 99% confidence 



level for one interesting parameter, xt 



6.63. We note for reference that a wavelength shift of 1 A corresponds to a velocity shift of 



62 km s 1 . The three examples shown here illustrate the following: Left: a significant shift of the offset peak of the broad line; Middle: 
subtle profile variation that mimics a shift of the broad line in the \ 2 curve; Right: no significant shift of the offset peak. 



TABLE 5 

Reasults of Cross-Correlation Analysis 
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The remaining objects have been omitted 



a Numbers with error bars denote a statistically significant shift. Pairs of numbers in square brackets 
denote upper and lower limits on the shift. A positive shift means that the post-SDSS spectrum is 
shifted towards longer wavelengths. 

b The broad H/3 profile varied substantially between observations, thus a meaningful shift could not 
be determined. These variaitons may include large changes in the profile shape or changes in the 
width that cause the two wings of the line to move in opposite directions. 

c Variability of the broad H/3 line profile affects our ability to determine a shift and increases the 
uncertainties. 

d The low S/N in the followup spectrum does not allow us to get meaningful constraints on the 
shift. 



variations in the broad H/3 profiles, such as small changes 
in the widths or variations in one of the two wings. Even 
though we do report shifts or limits in these cases, they 
are subject to an additional systematic error, which we 
have included in the limits we report in Table [5] In the 
remaining 16 cases we are not able to obtain meaningful 
measurements of velocity changes, or limits, either be- 
cause the broad H/3 profiles varied substantially between 
observations or because the S/N of of one of the two spec- 
tra was extremely low. We also note that iShields et al.l 
<|2009bf ) have determined an upper limit to the velocity 



change in J105041 of 31 ±60 km s" 1 (or [-30, 90], in the 
notation of Table [5]) . 

In Figure [10] we present our results graph ically (includ- 
ing th e limits on J105041, determined bv IShields et al.1 
I2009bl) . In the upper panel of this figure we plot the 
measured velocity change against the rest-frame time in- 
terval between observations. The statistically significant 
measurements are denoted by black, solid circles with 
99%-confidence error bars, while 99%-confidence limits 
on the velocity changes are shown as grey error bars. 
An alternative representation of the results is shown in 
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Fig. 10. — Top: Velocity change of broad H/3 peak (Ad) plotted 
against the rest-frame time interval between observations. Bottom: 
Accelerations determined from changes in the broad H/3 peak ve- 
locities (Au/At r est) plotted against the initial peak velocity (peak 
offset, as measured from the original SDSS spectra). In both pan- 
els, the black points with error bars show the objects where a sig- 
nificant change in the broad H/3 peak velocity was measured in 
the time interval covered by the observations. The light grey error 
bars show objects where only upper limits on the velocity change 
were obtained (inclu ded are also the limits derived for J105041 by 
IShields ct al. 20093)- All error bars and limits correspond to the 
99% (2.6ct) confidence level. 

the lower panel of the same figure where we plot the ob- 
served acceleration against the initial shift of the broad 
H/3 peak (as measured from the SDSS spectra; see Ta- 
ble SJ. We converted the velocity changes to accelera- 
tions (Ai>/Ai ros t, in km s _1 yr _1 ) by dividing the veloc- 
ity changes of Table [5] by the rest-frame time intervals 
between observations. The velocity changes appear to 
be distributed evenly with interval between observations. 
Moreover, there appears to be no relation between the 
acceleration and the initial velocity shift. We consider 
these findings further in our discussion of the results in 

m 

6.3. Limitations and Caveats 

As we noted above, there are a number of cases where 
the profile of the broad H/3 line has varied substantially 



between the SDSS and post-SDSS observations (see Ta- 
blets]). Such an effect represents an inevitable limitation 
to our ability to determine changes in the velocity of 
the broad lines, which relics on the long-term stability 
of the profiles. Variability of the line profiles can oc- 
cur on time scales comparable to the dynamical time of 
the BLR, which should be shorter than the orbital pe- 
riod of the hypothesized SBHB. Examples of objects with 
substantially variable broad H/3 profile include J112751, 
Jl 13706, J140700, and J180545 (see the comparison plots 
in Fig. [2]). A particularly interesting example of substan- 
tial variability is J093653, where the broad H/3 line de- 
clined in flux between the two epochs of observation with 
the result that the displaced peak of the broad line is now 
blended with the narrow line, making the measurement 
of a velocity change unreliable. 

More subtle profile variations can also affect our ability 
to determine velocity changes. For example, in J074157 
(see Fig. the broad H/3 line became narrower between 
the two observations and the line wings have moved 
closer to the line core, which prevents us from obtaining 
limits on the acceleration. In some cases, e.g., J091833, 
Jl 11916, and J154340 (see Fig. [2]), one of the two sides 
of the line has changed but the other has not; we can 
determine limits on the acceleration in such cases, but 
the uncertainties are larger than they would have been if 
the profile were stable, as we indicate in Table [5] 

Yet another type of profile variability that can influ- 
ence our results is a correlated change on the two sides 
of the broad H/3 profile (a seesaw-like pattern) that can 
mimic a velocity change. Such variability cycles have 
been observed in broad doub l e-pea ked emission lines 
(see, for example, iLewis et"al"1 12010). The beginnings 
of this pattern are discernible in the spectra of J031715 
and J160243 shown in Figure [21 for example. With this 
in mind, we note that some of the objects for which we 
report statistically significant velocity changes, specifi- 
cally those where the shifted peak of the broad H/3 line 
is partially blended with the narrow H/3 line, may be 
subject to this effect (e.g., J093844, J120924). Thus, the 
velocity changes reported in Table [5] are subject to ver- 
ification by future followup observations. There are, of 
course, additional reasons why further followup observa- 
tions are needed, as we discuss in more detail in the next 
section. 

7. DISCUSSION 

7.1. Alternative Interpretations of Displaced Emission 
Line Peaks 

Given that the starting point for our search for SB- 
HBs and rapidly-recoiling BHs was the selection of broad 
Balmer lines whose peaks are displaced from their nom- 
inal wavelengths, additional observations are needed to 
find the strong candidates from the initial list. The rea- 
son is that we cannot exclude alternative interpretations 
for the displaced peaks of the broad lines such as an un- 
usual structure or a perturbation of the BLR. For exam- 
ple, in the class of models where the BLR is the surface of 
the accretion disk or the base of the wind that is launched 
from it, perturbations in the form of large-scale spiral 
arms can create line profiles with a dominant peak or 
shoulder that is shifted from the nominal wavelength of 
the line. Examples of model profiles with such properties 
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can be found in lLewis et al.l (|2010t e.g., their Figure 36). 
In the same context, one can also envision cases where 
the broad Balmer lines from an unperturbed disk are 
double-peaked but a perturbation such as a spiral arm 
causes one of the two peaks to be considerably stronger 
than the other. Many examples of observed profiles that 
have such shapes c an be found in in the Appendix of 
IGezari et al.l (120071) : some of the more remarkable cases 
are reproduced in Figure 5 of lLauer fc Borosonl (|2009l) . 
Therefore, it is essential to carry out monitoring obser- 
vations to look for further variations in the velocities of 
the displaced peak of the broad line. 

The detection of velocity changes in broad H/3 lines 
between two epochs is only the hrst step in making the 
case that Balmer lines with displaced peaks are indeed 
signposts of SBHBs. Many monitoring observations are 
needed to verify that the velocity changes continue mono- 
tonically and that they follow the pattern expected for 
orbital motion. In particular the peaks of the broad 
line must drift from a negative velocity to a positive ve- 
locity and vice versa. Since our understanding of the 
structure and dynamics of the gas in the BLR is incom- 
plete, we cannot select SBHBs without observing com- 
plete orbital cycles because intrinsic variations of the 
broad line profiles on time scales of several years may 
resemble the pattern of orbital motion. A good case 
in point is provided by Mrk 668 (included in our sam- 
ple here as J140700). This object w as one of the first 
two proposed SBHBs ([Gaskelll 1 1983T ) and the variabil- 
ity of its broad Balmer lines has been monitored over 
a long period of time ( ~ 16 years; see IMarziani et al.l 
119931: IGezari et al.l [20071) . Although the velocity of the 
peak of the Balmer lines does change sign, the radial ve- 
locity curve is not compatible with orbital motion: the 
shape of the curve is not sinusoidal, the wings of the line 
pr ofile do not change accordingly (see Figures 30 and 31 
of IGezari et al.ll2007D . and the center of mass of a hy- 
pothesized binary ap pears to be at a diffe rent velocity 
than the host galaxy (|Marziani et al.l fl993T ) . This point 
is also bolstered by the case of 3C 390.3, a double-peaked 
emitter whose H/3 line changed in a manner remin iscent 
of orbital motion by ~ \ of a cycle (jGaskelllll996| ). The 
pattern changed drastically however aft er that and be- 
came incompatible with orbita l motion (|Eracleous et al.l 
Il997t iShapovalova et afll2001|) . More examples of broad 
Balmer lines with displaced peaks whose velocities have 
varied in a manner inconsistent wit h orbital mot i on ca n 
be found in IGezari et~aT1 (|2007| ) and lLewis et al.l (|2010f ). 
The need for long-term monitoring makes the search for 
SBHBs a challenging exercise in patience and persistence 
since the expected orbital periods can be quite long (see 
equation [4] and the associated discussion). 

Notwithstanding the caveats listed above, we consider 
the consequences of our findings for the two scenarios 
we have set out to test. In summary, the main obser- 
vational results of this work that bear directly on our 
tests are: (a) We have identified 88 quasars at z < 0.7 
(among 15,900 quasars from SDSS DR7) whose broad 
H/3 lines are shifted from their nominal positions by a 
fewx 1,000 km s _1 , (b) we have measured statistically 
significant accelerations of the offset peaks (in the range 
of -120 to +120 km s" 1 yr" 1 ) in 14/68 objects via fol- 
lowup observations, (c) we have measured the properties 



of the broad H/3 lines from the SDSS spectra and found 
a correlation between the skewness of the profiles and 
the velocity offset of the peak, which suggests that all 
the objects of this sample belong to the same family, i.e., 
have a single, common physical explanation and (d) we 
have measured the properties several narrow, forbidden, 
optical emission lines and found them to be very similar 
to those of typical quasars. 

7.2. Implications for Rapidly Recoiling BHs 

The rapidly recoiling BH hypothesis, remains a viable 
possibility for objects in this sample. This is because the 
primary observable, the velocity offset of the broad emis- 
sion lines, can be interpreted both as a result of an SBHB 
and as a result of a large recoil speed. Observational tests 
that have been proposed to distinguish between the two 
interpretations are indirect. The likelihood that a sig- 
nificant faction of our sample are rapidly recoiling black 
holes a ppears low a priori on theoretical grounds. In par- 
ticular, iDottreFIII (120101 ) find that only 0.2% of all recoil 
speeds should exceed 400 km s -1 , while the distribution 
of line shifts that we find, if all of these are interpreted 
as recoil speeds, implies that > 0.5% of recoil speeds ex- 
ceed 1,000 km s _1 . In the cases where we were able to 
measure an acceleration, a rapidly recoiling BH interpre- 
tation is not ruled out because intrinsic variability of the 
line profiles can mimic an acceleration as we discuss in 
N7.ll However, if the accelerations we have measured are 
confirmed to be real accelerations, that would weaken 
the case for a rapidly recoiling BH and strengthen the 
SBHB interpretation for those particular objects. 

An observational sele ction criterion f or rec oiling BH 
candidates, proposed by iBonning et al.l (|2007|) . is based 
on the widths of the narrow, forbidden lines. If the dis- 
placed quasar is illuminating the gas in its former host 
galaxy from the outside, the narrow-line region is not 
stratified in ionization, thus the widths of all the narrow, 
forbidden lines should be approximately the same. If, 
however, the gas emitting the narrow lines is associated 
with the recoiling BH (it could originate, for example, 
from bico nical outflows from th e accretion disk, as sug- 
gested bv lKomossa et all 120081 ). the widths of the lines 
would not be a good indicator of such a system since they 
may display all the terends observed in Seyfert galaxies. 
In such a case the luminosities and relative intensities of 
the narrow emission lines may provide a better test for 
recoiling BHs that their widths. 

Nevertheless, in the 50 quasars for which we were able 
to measure the FWHM of the [Ne v], [Ne ill], and [O ill] 
lines, there are 9 cases where the widths of these lines 
appear to be the same within 10%. In 8 of these cases, 
the [O n] lines are significantly narrower. The exception 
is J084313, where the width of the [O il] is compara- 
ble to that of the other forbidden lines. In this object, 
the velocity offset of the broad H/3 line from the system 
of the narrow lines is only —610 km s _1 and the pro- 
file of the line if fairly symmetric, as indicated by the 
skewness coefficients (this is another important consid- 
eration for selecting candidate recoiling BHs according to 
IShields et al.l 12009a ). In spite of the above, the relative 
strengths of the narrow, forbidden lines of this object 
are very typical of a Seyfert galaxy, which suggests that 
the ionization parameter of the gas in the narrow-line 
region should be similar to that of Seyfcrts. This is not 
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in agreement with the proposed picture of a displaced 
quasar that illuminates its host galaxy from a significant 
distance. The quasar J0 92712, proposed as a candidate 
rapidly recoiling BH by iKomossa et all ()2008f ). has two 
systems of narrow lines and in the one that would be 
attributed to the (former) host galaxy in this interpre- 
tation the narrow lines are unusually narrow and have 
similar widths. How ever, alternative interpretations have 
also been proposed (iBogdanovic et al . 2009; Dotti et "all 
I2009at [Shields et al.ll2009at iHeckman et al.ll2009t ). which 
means that the rapidly recoiling BH scenario is not a 
unique interpretation for this very unusual object. 

We note in conclusion that the process of photoioniza- 
tion of the gas in a galaxy by an offset or neighboring 
AGN warrants further investigation in order to develop 
more quantitative and robust emission-line diagnostics. 
O n the th e oretic al side, the photoionization calculations 
of lGnedinl (fl997h can be extended to cover a wider range 
of conditions and to make predictions for the relative 
intensities, luminosities, and profiles of the narrow emis- 
sion lines. On the observational side, progress can be 
made by seeki ng out and studyin g more analogs of the 
Was 49 system (jMoran et al.lll992l ) , a dual Seyfert galaxy 
in which the more luminous AGN contributes to the ion- 
ization of the interstellar medium of the host galaxy of 
its less luminous neighbor. 

7.3. Implications for SBHBs 

Considering our results in the context of the SBHB 
hypothesis, we ca n compare them with the specific pre- 
dictions made by IVolonteri et al.l (|2009() for the number 
of such systems that may be present in the SDSS DR7 
quasar sample at z < 0.7. If we suppose that all of the 
objects in our sample are SBHBs, then their number is 
broadly consistent with one of the two models consid- 
ered bv IVolonteri et afl <|2009f ) . their model I. According 
to this model, which is based on the assumption that 
all SBHBs are active (i.e., accreting) regardless of the 
merger history of their host, there are at most 160 sys- 
tems with q > 10~ 2 in the SDSS DR7 quasar sample. 
If we only take the 14 systems with detectable velocity 
changes to be SBHBs, their number i s broa dly consis- 
tent with model II of IVolonteri et al.l (|2009() according 
to which SBHBs are active only in the aftermath of a 
major merger. Under the assumptions of this model 
the predicted number of SBHBs with q > 10~ 2 in the 
SDSS DR7 quasar sample is 16. The comparison of the 
above predictions with the observational results is com- 
plicated by a number of uncertainties. On the theoretical 
side: (a) the models refer to a volume-limited sample of 
quasars, while the SDSS quasar sample is not volume lim- 
ited (given the flux limit of the quasar catalog, quasars 
with i magnitudes between 19.1 and 21.2 are missed even 
though they would sattisfy the quasar luminosity crite- 
rion), (b) the predicted numbers quoted above rely on 
the assumption that both of the BHs in a binary may 
accrete (with an assumed distribution of Eddington ra- 
tios), without regard to their relative luminosity, while 
our selection method picks out systems in which the BLR 
around only one of the two BHs is visible, (c) inclusion 
in the theoretical sample is based on the luminosity of 
the primary, regardless of whether or not it is less lumi- 
nous than the secondary. The above uncertainties can 
change the upper limits quoted above in cither direction. 



On the observational side: (a) some of the candidates we 
report may turn out not to be SBHBs for the reasons 
discussed in H7.1[ (b) if we assume that our selection 
method is appropriate for finding SBHBs, we are likely 
to have missed systems with small projected orbital ve- 
locities (i.e., systems viewed close to face on, or caught 
near conjunction, or with a combination of low masses 
and wide separations), and (c) the number of objects 
for which a measurable acceleration is likely to change 
change with future observations. 

Assuming that velocity changes we have measured cor- 
respond to orbital accelerations in an SBHB we can ex- 
plore the implications of the results for this hypothesis. 
The correlation between the skewness of the broad H/3 
profiles and the shift of their peaks (see Figure [5]) sug- 
gests that all the objects in this sample belong to the 
same familjQ, therefore we ask why an acceleration was 
measurable only in a small fraction of the sample. As 
we noted in 3Hwe expect that the SBHBs we can detect 
with our strategy can have periods from a few decades 
to a few centuries depending on the total mass of the 
binary and the observed velocity shift, which constrains 
the orbital separa tion given th e masses of the two BHs. 
As pointed out bv lLoebl (|2010| ). the likelihood of finding 
a system with a given set of intrinsic properties depends 
on its decay time via emission of gravitati onal radiation . 
This can be estimated from the formula of lPetersI (|1964l) , 
which we can cast in terms of intrinsic system properties 
and observables using equation (j4j as follows: 
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Since the lifetime of an observed SBHB is such a sensi- 
tive function of the observed velocity offset of the broad 
emission line, 112,3, wc would expect that low- mass, high- 
speed binaries should be extremely rare since their or- 
bits should be rapidly decaying by gravitational radia- 
tion. Specifically, setting Mg, = 0.1 and 1*2,3 — 2 in 
equation ([7j gives t gI «3x 10 8 yr, while setting Mg = 1 
and U2.3 = 3 gives i gr w 1 X 10 s yr. The above consid- 
erations lead us to conclude that systems with observed 
velocity shifts in excess of 1,000 km s _1 may be (a) highly 
inclined (hence with large separations, according to equa- 
tion [J), (b) caught near quadrature, (c) very massive, 
or (d) very small mass-ratio systems (q <C 0.1). Under 
different conditions, such systems would have very short 
lifetimes. The relatively uniform distribution of velocity 
changes and accelerations depicted in Figure ITUl suggests 
that a combination of the above four effects is at work. 
Highly inclined or massive systems would have low accel- 
erations, according to equation ([5]), which may explain 
why we were not able to detect an acceleration in the 
majority of the objects we have followed up. Continued 

12 This is an optimistic hypothesis. A more pessimistic interpre- 
tation of this correlation is, of course, that none of the objects in 
our sample are SBHBs. 
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monitoring can discriminate among the above possibili- 
ties by providing constraints on the binary phase through 
the shape of the radial velocity curve. 

Useful constraints on the properties of a hypothesized 
SBHB can be obtained from a few observations span- 
ning a baseline of about a decade. Accurate measure- 
ments of the velocities of the peaks over such a base- 
line can constrain the curvature of the radial velocity 
curve and yield a lower limit on the period and to- 
tal mass of the SBHB. Th is technique, pioneered by 
Halpern fc Filippenkol (|1988j ) and developed further by 
Eracleous et al.l (|1997l ). was used to test and eventu- 
ally reject the SBHB hypothesis for three quasars with 
double-peaked Balmer lines. An application of the same 
technique to objects from the present sample can lead 
to similar constraints, which could render the SBHB 
hypothesis untenable if these constraints are restrictive 
enough. To this end, we are continuing to monitor all of 
our targets with the immediate goal of collecting a few 
radial velocity measurements for each object and using 
them to evaluate the plausibility of the SBHB hypothesis 
by the above method. 

Another avenue for testing the SBHB hypothesis for 
shifted broad Balmer line peaks involves a comparison 
of the profiles of the Lya and Balmer lines. If the 
displaced peaks originate in a perturbed disk around a 
single BH rather than in a BLR bound to one of the 
two BHs in SBHB, we would expect the profiles of the 
broad optical and UV emission lines to differ consid- 
erably This expectation is based on an analogy with 
nearby, well-studied double-peaked emitters (see, for ex- 
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